Two-Qubit State Tomography using a Joint Dispersive Read-Out 
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Quantum state tomography is an important tool in quantum information science for complete 
characterization of multi-qubit states and their correlations. Here we report a method to perform a 
joint simultaneous read-out of two superconducting qubits dispersively coupled to the same mode of 
a microwave transmission line resonator. The non-linear dependence of the resonator transmission 
on the qubit state dependent cavity frequency allows us to extract the full two-qubit correlations 
without the need for single shot read-out of individual qubits. We employ standard tomographic 
techniques to reconstruct the density matrix of two-qubit quantum states. 
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Quantum state tomography allows for the reconstruc- 
tion of an a-priori unknown state of a quantum system 
by measuring a complete set of observables [1]. It is an 
essential tool in the development of quantum informa- 
tion processing [2] and has first been used to reconstruct 
the Wigner-function [3[ of a light mode [4] by homodyne 
measurements, as suggested in a seminal paper by Vogel 
and Risken [5] . Subsequently, state tomography has been 
applied to other systems with a continuous spectrum, 
for instance, to determine vibrational states of molecules 
0], ions and atoms [8]. Later, techniques have been 
adapted to systems with a discrete spectrum, for exam- 
ple nuclear spins Q , polarization entangled photon pairs 
[lOj], electronic states of trapped ions [11], states of hy- 
brid atom-photon systems [l2[, and spin-path entangled 
single neutrons [131 ]. 

Recent advances have enabled the coherent control 
of individual quantum two-level systems embedded in 
a solid-state environment. Numerous experiments have 
been performed with superconducting quantum devices 
manifesting the rapid progress and the promising 
future of this approach to quantum information process- 
ing. In particular, the strong coupling of superconducting 
qubits to a coplanar waveguide resonator can be exploited 
to perform cavity quantum electrodynamics (QED) ex- 
periments on a chip [H, [H, [l7| in an architecture known 
as circuit QED. The high level of control over the dy- 
namics of this coupled quantum system has been demon- 



strated, e. g., in [18|, [19[. State tomographic methods 



have already been used in superconducting circuits to 
verify the entanglement between two phase qubits [20j . 
There, the state is determined for each individual qubit 
with single-shot read-out such that two-qubit correlations 
can be evaluated by correlating the single measurement 
outcomes. In contrast, in this letter we extract two- 
qubit correlations from the simultaneous averaged mea- 
surement of two qubits dispersively coupled to a common 
resonator. This possibility has also been pointed out in 



Ref.[21. 

In the setup shown in Fig. [IJ two superconducting 
qubits are coupled to a transmission line resonator op- 
erating in the microwave regime [ll|. Due to the large 
dipole moment of the qubits and the large vacuum 
field of the resonator the strong coupling regime with 
<7i,2 ^> ft>7i is reached. #i/27r « #2/27r = 133 MHz de- 
notes the similar coupling strenghts of both qubits and 
k/2tt w 1.65 MHz, 7i/2tt w 0.25 MHz the photon and the 
qubit decay rates, respectively. The qubits are realized 
as transmons [22], a variant of a split Cooper pair box 
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FIG. 1: Schematic of the experimental setup with two qubits 
coupled via the capacitances C g to a microwave resonator op- 
erated at a temperature of about 20 mK. The transition fre- 
quencies of the qubits are adjusted by external fluxes $i and 
$2. The resonator-qubit system is probed through the input 
and output capacitances Cm and C ou t by a microwave sig- 
nal at frequency ou m • Additionally, local control of the qubits 
is implemented by capacitively coupled signals Udi and cdd2, 
which are phase and amplitude modulated using IQ mixers. 
The output signal is detected in a homodyne measurement at 
room temperature. 
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|23( | with exponentially suppressed sensitivy to 1 /f charge 
noise [24j. The transition frequencies uo a j (j = 1, 2) of the 
qubits are tuned separately by external flux bias coils. 
Both qubits can be addressed individually through lo- 
cal gate-lines using amplitude and phase modulated mi- 
crowaves at frequencies uodi and cod2- Read-out is ac- 
complished by measuring the transmission of microwaves 
applied to the resonator input at frequency oj m close to 
the fundamental resonator mode uo r . At large detunings 
Aj = uu a j —uo r of both qubits from the resonator, the dis- 
persive qubit-resonator interaction gives rise to a qubit 
state dependent shift of the resonator frequency. In this 
dispersive limit and in a frame rotating at cj m the rele- 
vant Hamiltonian reads 251 



') -(x,+yj (x-x z ) (tl 2 -x,) 



H = h(A rrn + xiVzi+X2Vz2)a t a (1) 
I" Xj)<?zj + he{t){a ] + a), 



where A rm = a; r — a; m is the detuning of the measure- 
ment drive from the resonator frequency. The coefficients 
Xi,2 are determined by the detuning Ai ? 2, the coupling 
strength g\^ and the design parameters of the qubit [221 ]. 
The last term in Eq. (pQ) models the measurement drive 
with amplitude e(t). 

The operator x = Xi^i + X2&z2, which describes 
the dispersive shift of the resonator frequency, is lin- 
ear in both qubit states. It does not contain two- 
qubit terms like cr z \cr Z 2 from which information about 
the qubit-qubit correlations could be obtained. How- 
ever, in circuit QED instead of measuring frequency 
shifts directly, we record quadrature amplitudes of mi- 
crowave transmission through the resonator which de- 
pend nonlinear ly on these shifts. The average values 
of the field quadratures (/(£)) = Tr [/)(£) (at + a)] and 
(Q) = iTr\p(t)(aJ — a)] are determined from the am- 
plified voltage signal at the resonator output in a ho- 
modyne measurement. Here, p(t) = Um(t)p(0)Um(ty 
denotes the time evolved state of both resonator and 
qubit under measurement. In the dispersive approxima- 
tion we can safely assume this state to be separable be- 
fore the measurement, which is taken to start at time t m , 
p(tm) — |0)(0| &pq(t m ). Using these expressions, we find 
(Q(t)) =iTV,[ / 5 q (i ro )(0|C/t i (i)(at-a)U- TO (f)|0)] (and sim- 
ilarly for (J(t))), where Tr q denotes the partial trace over 
the qubit. This expression is evaluated using the input- 
ouput formalism [26[ including cavity decay k. In the 
steady-state, this yields (J) s , (Q) s = -e Tr g [^(^M^q] 
with 



Mi = 
M Q = 



2(A rm + x) 



(A rm + x) 2 + (^/2) 2 



(A rm + x) 2 + (^/2) 



2 ' 



(2) 
(3) 



We note that the measurement operators are nonlinear 
functions of \. Thus, M/ ? q comprises in general also two- 
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FIG. 2: (a) Q-quadrature of the resonator field for the qubits 
in states gg, eg, ge and ee as a function of the detuning 
A r m- Tomography measurements have been performed at 
A rm = (xi + X2) indicated by an arrow, (b) Measured (data 
points) time evolution of the Q-quadrature for the indicated 
initial states compared to numerically calculated responses 
(solid lines). All parameters have been determined in inde- 
pendent measurements. 



qubit correlation terms proportional to (J z \Vz2, which al- 
low to reconstruct the full two-qubit state. 

In our experiments the phase of the measurement mi- 
crowave at frequency A rm = (xi + X2) is adjusted such 
that the Q-quadrature of the transmitted signal carries 
most of the signal when both qubits are in the ground 
state. The corresponding measurement operator can be 
expressed as 



M 



id + (3 w a z i + (3 ia z2 + j3na zl a z2 



(4) 



where — a__ + (- 
with the coefficients 



-l)>a.+ + (-!)*<*+. + {-l) i+ ia ++ , 



a±± = -€k{(k/2) 2 + (Arm ± Xi ± X2) 2 } 



-1/2 



(5) 



representing the qubit state dependent Q-quadrature am- 
plitudes of the resonator field in the steady-state limit 
and for infinite qubit-lifetime (Fig. [2fa)). 

Since we operate in a regime, where the qubit relax- 
ation cannot be neglected, the steady-state expression 
is of limited practical use. The decay of a qubit to 
its ground state changes the resonance frequency of the 
resonator and consequently limits the read-out time to 
~ I/71. A typical averaged time-trace of the resonator 
response for pulsed measurements is shown in Fig. Efb), 
similar to the data presented in [21]. The qubits are 
prepared initially in the states |ee), \eg), \ge) and \gg), 
respectively, using the local gate lines. The time depen- 
dence of the measurement signal is determined by the rise 
time of the resonator and the decay time of the qubits. 
It is in excellent agreement with calculations (solid lines 
in Fig. 12(b)) °f tne dynamics of the dispersive Jaynes- 



Cummings Hamiltonian [25|, [2 



Due to the quantum non-demolition nature of the 
measurement [l6[, M remains diagonal in the instan- 
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taneous qubit eigenbasis during the measurement pro- 
cess. Therefore, a suitable realistic measurement opera- 
tor M' can be defined by replacing the a±± in Eq. ([5]) 
with the integrated signal from t m to the final time T, 
a'±± = 1/N ((M(t)) ±± -(M(t))..)dt with the ground 

state response (M (£))__ subtracted. The normalization 
constant N is chosen such that a' + _ = 1. 

To reconstruct the combined state p q of both qubits, 
a suitable set of measurements has to be found to deter- 
mine unambiguously the 16 coefficients nj of the density 
matrix p q — Ylij=o r ij &i ® &j with the identity &q = id 
and <72, (J3} = {a x , a y , a z }. Such a complete set 
of measurements is constructed by applying appropriate 
single qubit rotations G SU(2) SU(2) before the 
measurement in order to measure the expectation val- 
ues (M k ) = Ti[MU k p q Ui] = Tr[ulMU k p q ]. The latter 
equality defines the set of measurement operators = 
UlMUk- This illustrates again that a measurement op- 
erator M involving non-trivial two-qubit terms (Tixcr^ is 
necessary for state tomography. In fact, single-qubit op- 
erations Uk = Uki<S>Uk2 alone cannot be used to generate 
correlation terms since Ul (id §§a z )Uk = id0 (Ul 2 a z Uk2), 
for instance. As Tr [(cr^ <g> cri)(cr m <g> a n )] = Sk m Si n , some 
coefficients of the density matrix p q would not be de- 
termined in an averaged measurement. 

To identify the coefficients we perform 16 linearly- 
independent measurements. The condition for the com- 
pleteness of the set of tomographic measurements is the 
non-singularity of the matrix A defined by the relation 
(Mk) = J2ilo^klTi between the the expectation values 
(Mk) and the coefficients of the density matrix r\ with 
I = i + 4j. This condition is only violated if one of the 
coefficients of M in Eq. (j4]) vanishes. For instance, 
An = Alo = for A rm = 0, which reflects the fact that 
we cannot distinguish two identical qubits due to sym- 
metry reasons as apparent from Fig. E^a). 

Our pulse scheme for the state tomography is shown 
in Fig. [3l The transition frequencies of the qubits are 
adjusted to o; a i/27r = 4.5 GHz and co a 2/27T = 4.85 GHz. 
At this detuning from the resonator frequency o; r /27r = 
6.442 GHz the cavity pulls are xi — ~ 1 MHz and 
X2 = —1.5 MHz [22]. First, a given two-qubit state 
is prepared. Then a complete set of tomography mea- 
surements is formed by applying the combination of 
{(71-/2)3., (7r/2) y , (7r),id} pulses to both qubits over their 
individual gate lines using amplitude and phase con- 
trolled microwave signals. The wanted rotation angles 
are realized with an accuracy better than 4°. Finally, 
the measurement drive is applied at uo m = 6.445 GHz 
corresponding to the maximum transmission frequency 
of the resonator with both qubits in the ground state. 
To determine the measurement operator M 7 , 7r-pulses 
are alternately applied to both qubits to yield signals 
as shown in Fig. [21(b). From this data the coefficients 
(P' 00 , (3' 01 ,(3' w , (3' n ) = (0.8, -0.3, -0.4, -0.1) of M' are de- 
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FIG. 3: Pulse scheme for state tomography, see text. 



duced, where the non- vanishing /3' u allows for a mea- 
surement of arbitrary quantum states. As an exam- 
ple of this state reconstruction, in Fig. ID[a) the ex- 
tracted density matrix p q of the product state |^ S ep) = 
1/V2(\g) + \e))®l/V2(\g) + i|e» is shown. In Fig.^b) 
the Bell state |$) = l/y/2(\g) \g) - i\e) <%> \e)) prepared 
by a sequence of sideband pulses [25|, 0, [23] is recon- 
structed. 6.6 x 10 4 and 6.6 x 10 5 records have been 
averaged, respectively, for each of the 16 tomographic 
measurement pulses to determine the expectation val- 
ues (M' k ) for the two states. The corresponding ideal 
state tomograms are depicted in Fig. \Ms) and (d). To 
avoid unphysical, non positive-semidefinite, density ma- 
trices originating from statistical uncertainties, all to- 
mography data has been processed by a maximum like- 
lihood method [13, The corresponding fidelities 
= (Wp q W) 1/2 are ^ sep = 95% and T*- = 74%. 
These results are in close agreement with theoretically 
expected fidelities when taking finite photon and qubit- 
lifetimes into account. As a result, the loss in fidelity and 
concurrence of the Bell state are not due to measurement 
errors but to the long preparation sequence [29[ . 

In conclusion, we have presented a method to jointly 
and simultaneously read-out the full quantum state of 
two qubits dispersively coupled to a microwave resonator. 
In a measurement of the field quadrature amplitudes of 
microwaves transmitted through the resonator each pho- 
ton carries information about the state of both qubits. In 
this way the two-qubit correlations can be extracted from 
an averaged measurement of the transmission amplitude 
without the need for single shot or single qubit read-out. 
This method can also be extended to multi-qubit systems 
coupled to the same resonator mode. 
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FIG. 4: Real and imaginary part of reconstruced density 
l/y/2(\g) + i|e)) and (b) the Bell state |$) = l/v^2 (|p> <g> |p) 
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